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1. NEUTRINO MASS, MIXING, 
AND FLAVOR CHANGE* 

Revised March 2008 by B. Kaysert 

Fermilab, P.O. Box 500, Batavia, IL 60510, USA 

There is now compelling evidence that atmospheric, solar, accelerator, and reactor 
neutrinos change from one flavor to another. This implies that neutrinos have masses and 
that leptons mix. In this review, we discuss the physics of flavor change and the evidence 
for it, summarize what has been learned so far about neutrino masses and leptonic 
mixing, consider the relation between neutrinos and their antiparticles, and discuss the 
CN ■ open questions about neutrinos to be answered by future experiments. 

/. The physics of flavor change: If neutrinos have masses, then there is a spectrum 
of three or more neutrino mass eigenstates, v\, z^2, ^3, • • -, that are the analogues of 
the charged-lepton mass eigenstates, e, fj,, and r. If leptons mix, the weak interaction 
coupling the W boson to a charged lepton and a neutrino can couple any charged-lepton 
mass eigenstate £ a to any neutrino mass eigenstate i/j. Here, a = e, //, or r, and £ e is the 
electron, e£c.The amplitude for the decay of a real or virtual W + to yield the specific 
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combination £+ + v^ is U*^, where U is the unitary leptonic mixing matrix [1]. Thus, the 
<-| \ neutrino state created in the decay W + — > £^ + v is the state 

^. Wa)=Y, U aM) ■ (1-1) 

>, 

This superposition of neutrino mass eigenstates, produced in association with the charged 
lepton of "flavor" a, is the state we refer to as the neutrino of flavor a. Assuming CPT 
invariance, the unitarity of U guarantees that the only charged lepton a v a can create in a 
detector is an £ a , with the same flavor as the neutrino. Eq. (1.1) may be inverted to give 
O ■ v^ 

Wi) = Z^UpiM , (1-2) 


which expresses the mass eigenstate Vi as a superposition of the neutrinos of definite 
^% ' flavor. 

While there are only three (known) charged lepton mass eigenstates, it may be that 
there are more than three neutrino mass eigenstates. If, for example, there are four z^, 
then one linear combination of them, 

\"») =T,Uri\"i) > (1-3) 

i 

does not have a charged-lepton partner, and consequently does not couple to the Standard 
Model W boson. Indeed, since the decays Z — > v a V a of the Standard Model Z boson 
have been found to yield only three distinct neutrinos v a of definite flavor [2], v s does 
not couple to the Z boson either. Such a neutrino, which does not have any Standard 
Model weak couplings, is referred to as a "sterile" neutrino. 

Neutrino flavor change is the process v a — > ug, in which a neutrino born with flavor 
a becomes one of a different flavor (3 while propagating in vacuum or in matter. This 
process, often referred to as neutrino oscillation, is quantum mechanical to its core. 
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2 1 . Neutrino mixing 

Rather than present a full wave packet treatment [3], we shall give a simpler description 
that captures all the essential physics. We begin with oscillation in vacuum, and work in 
the neutrino mass eigenstate basis. Then the neutrino that travels from the source to the 
detector is one or another of the mass eigenstates v^. The amplitude for the oscillation 
v a - ► v 3i Amp [y a — ► vg), is a coherent sum over the contributions of all the Vi, given by 

Amp (u a -> vp) = Y, Ki Pr °P iy%) U (3i • ( 1A ) 



In the contribution U*^ Prop (i/j) Uqi of vi to this sum, the factor U^ is the amplitude 
for the neutrino v a to be the mass eigenstate z/j [see Eq. (1.1)], the factor Prop [yi) is 
the amplitude for this z/j to propagate from the source to the detector, and the factor 
Uqi is the amplitude for the v^ to be a vq [see Eq. (1.2)]. From elementary quantum 
mechanics, the propagation amplitude Prop [yi) is exp[— imjTj], where m$ is the mass of 
vi, and Ti is the proper time that elapses in the Vi rest frame during its propagation. 
By Lorentz invariance, rrnTi = E{t — p%L, where L is the lab-frame distance between the 
neutrino source and the detector, t is the lab-frame time taken for the beam to traverse 
this distance, and E^ and pj are, respectively, the lab-frame energy and momentum of the 
Vi component of the beam. 

In the probability P(v a — ► vg) = | Amp (y a — >■ vg)\ for the oscillation v a — >• z/^, only 
the relative phases of the propagation amplitudes Prop (z^) for different mass eigenstates 
will have physical consequences. From the discussion above, the relative phase of Prop (z/j) 
and Prop (vj), Sfaj, is given by 

p i -p j )L-(E i -E j )t . (1.5) 



In practice, experiments do not measure the transit time t. However, Lipkin has shown [4] 
that, to an excellent approximation, the t in Eq. (1.5) may be taken to be L/v, where 

Pi+Pj n ^ 

Ei + Ej 

is an approximation to the average of the velocities of the z/j and Vj components of the 
beam. Then 



r^j 



\J\Uf)n J-/ J_/ I I IV A I I VA I „ ^_ 

J Pi+Pj Pi+Pj V J */ 2£ 

where, in the last step, we have used the fact that for highly relativistic neutrinos, 
Pi and pj are both approximately equal to the beam energy E. We conclude that 
all the relative phases in Amp (v a — > vq\ Eq. (1.4), will be correct if we take 
Prop (z/j) = exp (—im?L/2E), so that 

Amp („ a - I/,,) = Y, Ki e~^ L ' 2E U Bl . (1.8) 
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1. Neutrino mixing 3 

Squaring, and making judicious use of the unitarity of U, we then find that 



p( 



Vo 



■ Vp) = $a[3 

4£» (KiUpiUajU^j) sin 2 [l.27Am| (L/E) 



i>j 



+2^3 [K t U pi U aj U^ sin [2.54 Am% (L/E) 
i>j 



(1.9) 



o 

Here, Am^ = m^ — m| is in eV , L is in km, and E is in GeV. We have used the fact 
that when the previously omitted factors of h and c are included, 



Amfj (L/AE) ~ 1.27 Am^- (eV' 



L(km) 

£ (GeV) 



Assuming that CPT invariance holds, 

P (v a ->vp) = P {yp 

But, from Eq. (1.9) we see that 

P{v (3 ^v a -U) = P{v a - 
Thus, when CPT holds, 

PiVa^Vp-U) =P[y a - 



'0 
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(1.10) 



(i.ii) 



(1.12) 



(1.13) 



That is, the probability for oscillation of an antineutrino is the same as that for a 
neutrino, except that the mixing matrix U is replaced by its complex conjugate. Thus, if 
U is not real, the neutrino and antineutrino oscillation probabilities can differ by having 
opposite values of the last term in Eq. (1.9). When CPT holds, any difference between 
these probabilities indicates a violation of CP invariance. 

As we shall see, the squared-mass splittings Am 2 - called for by the various reported 
signals of oscillation are quite different from one another. It may be that one splitting, 
AM , is much bigger than all the others. If that is the case, then for an oscillation 
experiment with L/E such that AM L/E = 0(1), Eq. (1.9) simplifies considerably, 
becoming 

1.27 AM 2 (L/E) 



PlV a ^%) ^S a(3 sin 2 



for (3 7^ a, and 



PH ~^ 



is a )~l-4T a (l-T a )sm 2 1.27 AM 2 (L/E) 



Here, 
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a/3 



i Up 



(1.14) 



(1.15) 



(1.16) 
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4 1 . Neutrino mixing 

and 



T a =J2 Pail 2 > (1-17) 

i Up 



where "i Up" denotes a sum over only those neutrino mass eigenstates that lie above 
AM or, alternatively, only those that lie below it. The unitarity of U guarantees that 
summing over either of these two clusters will yield the same results for S a p and for 
T a (l-T a ). 

The situation described by Eqs. (1.14)-(1.17) may be called "quasi-two-neutrino 
oscillation." It has also been called "one mass scale dominance" [5]. It corresponds to 
an experiment whose L/E is such that the experiment can "see" only the big splitting 
AM . To this experiment, all the neutrinos above AM appear to be a single neutrino, 
as do all those below AM . 

The relations of Eqs. (1.14)-(1.17) apply to a three-neutrino spectrum in which one of 
the two squared-mass splittings is much bigger than the other one. If we denote by z/3 the 
neutrino that is by itself at one end of the large splitting AM , then S a p = 4\U a sU^\ 
and T a = \U a s\ . Thus, oscillation experiments with AM L/E = 0(1) can determine 
the flavor fractions |t/ a 3| of z/3. 

The relations of Eqs. (1.14)-(1.17) also apply to the special case where, to a good 
approximation, only two mass eigenstates, and two corresponding flavor eigenstates (or 
two linear combinations of flavor eigenstates), are relevant. One encounters this case 
when, for example, only two mass eigenstates couple significantly to the charged lepton 
with which the neutrino being studied is produced. When only two mass eigenstates 
count, there is only a single splitting, Am , and, omitting irrelevant phase factors, the 
unitary mixing matrix U takes the form 



U 





v\ 


"2 


V a 


cos 9 


sin0 


V(3 


— sin 9 


cos 9 



(1.18) 



Here, the symbols above and to the left of the matrix label the columns and rows, and 9 is 
referred to as the mixing angle. From Eqs. (1.16) and (1.17), we now have S a p = sin 29 
and 4T a (l — T a ) = sin 20, so that Eqs. (1.14) and (1.15) become, respectively, 

P ( ( P a ^vp) = sin 2 29 sin 2 1.27 Am 2 (L/E) (1.19) 

with (3 7^ a, and 

P ( c ^ -> ( P tt ) = 1 - sin 2 29 sin 2 [l .27 Am 2 (L/E)] . (1.20) 

Many experiments have been analyzed using these two expressions. Some of these 
experiments actually have been concerned with quasi-two-neutrino oscillation, rather than 
a genuinely two-neutrino situation. For these experiments, "sin 20" and "Am have the 
significance that follows from Eqs. (1.14)-(1.17). 
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1. Neutrino mixing 5 

When neutrinos travel through matter (e.g., in the Sun, Earth, or a supernova), their 
coherent forward-scattering from particles they encounter along the way can significantly 
modify their propagation [6]. As a result, the probability for changing flavor can be 
rather different than it is in vacuum [7]. Flavor change that occurs in matter, and that 
grows out of the interplay between flavor-nonchanging neutrino-matter interactions and 
neutrino mass and mixing, is known as the Mikheyev-Smirnov-Wolfenstein (MSW) effect. 

To a good approximation, one can describe neutrino propagation through matter via a 
Schrodinger-like equation. This equation governs the evolution of a neutrino state vector 
with several components, one for each flavor. The effective Hamiltonian in the equation, a 
matrix Ti in neutrino flavor space, differs from its vacuum counterpart by the addition of 
interaction energies arising from the coherent forward neutrino-scattering. For example, 
the u e -u e element of Ti includes the interaction energy 

V = V2G F N e , (1.21) 

arising from VF-exchange-induced v e forward-scattering from ambient electrons. Here, Gp 
is the Fermi constant, and N e is the number of electrons per unit volume. In addition, the 
u e —u e , v^—ffx, and v T -v T elements of Ti all contain a common interaction energy growing 
out of Z-exchange-induced forward-scattering. However, when one is not considering the 
possibility of transitions to sterile neutrino flavors, this common interaction energy merely 
adds to Ti a multiple of the identity matrix, and such an addition has no effect on flavor 
transitions. 

The effect of matter is illustrated by the propagation of solar neutrinos through 
solar matter. When combined with information on atmospheric neutrino oscillation, the 
experimental bounds on short-distance (L ^ 1 km) oscillation of reactor V e [8] tell us 
that, if there are no sterile neutrinos, then only two neutrino mass eigenstates, v\ and ^2, 
are significantly involved in the evolution of the solar neutrinos. Correspondingly, only 
two flavors are involved: the v e flavor with which every solar neutrino is born, and the 
effective flavor v x -- some linear combination of v^ and v T — which it may become. The 
Hamiltonian Ti is then a 2 x 2 matrix in u e —u x space. Apart from an irrelevant multiple 
of the identity, for a distance r from the center of the Sun, Ti is given by 



H =Ti v + Ti M (r) 



ArriQ 
AE 



— cos 26*0 sin 26*0 
sin 26*0 cos 26*© 



+ 



V(r) 




;i.22) 



Here, the first matrix Tiy is the Hamiltonian in vacuum, and the second matrix TiM(r) 
is the modification due to matter. In Tiy, 9@ is the solar mixing angle defined by the 
two-neutrino mixing matrix of Eq. (1.18) with 9 = 9@, v a = u e , and up = v x . The 
splitting Airig is m^ — Tn\, and for the present purpose we define u<i to be the heavier 
of the two mass eigenstates, so that Am is positive. In TiM( r ), V(r) is the interaction 
energy of Eq. (1.21) with the electron density N e (r) evaluated at distance r from the 
Sun's center. 

From Eqs. (1.19-1.20) (with 9 = 9q), we see that two-neutrino oscillation in vacuum 
cannot distinguish between a mixing angle 9q and an angle 9' & = 7r/2 — 9q. But these 
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6 1 . Neutrino mixing 

two mixing angles represent physically different situations. Suppose, for example, that 
#0 < 7r/4. Then, from Eq. (1.18) we see that if the mixing angle is #©, the lighter 
mass eigenstate (defined to be u\) is more v e than u x , while if it is #q, then this 
mass eigenstate is more v x than v e . While oscillation in vacuum cannot discriminate 
between these two possibilities, neutrino propagation through solar matter can do so. 
The neutrino interaction energy V of Eq. (1.21) is of definite, positive sign [9]. Thus, the 
v e -v e element of the solar 7i, — (AitIq/AE) cos 20© + V(r) : has a different size when the 
mixing angle is 0'q = n/2 — 9@ than it does when this angle is #0. As a result, the flavor 
content of the neutrinos coming from the Sun can be different in the two cases [10]. 

Solar and long-baseline reactor neutrino data establish that the behavior of solar 
neutrinos is governed by a Large-Mixing- Angle (LMA) MSW effect (see Sec. II). Let us 
estimate the probability P{y e —> v e) that a solar neutrino that undergoes the LMA-MSW 
effect in the Sun still has its original v e flavor when it arrives at the Earth. We focus 
on the neutrinos produced by B decay, which are at the high-energy end of the solar 
neutrino spectrum. At r ~ 0, where the solar neutrinos are created, the electron density 
N e ~ 6 x 10 /cm [11] yields for the interaction energy V of Eq. (1.21) the value 
0.75 x 10 -5 eV 2 /MeV. Thus, for Am 2 ^ in the favored region, around 8 x 10 -5 eV 2 , 
and E a typical B neutrino energy (~6-7 MeV), Hm dominates over Tiy. This means 
that, in first approximation, TC(r ~ 0) is diagonal. Thus, a B neutrino is born not 
only in a v e flavor eigenstate, but also, again in first approximation, in an eigenstate 
of the Hamiltonian 7i(r ~ 0). Since V > 0, the neutrino will be in the heavier of the 
two eigenstates. Now, under the conditions where the LMA-MSW effect occurs, the 
propagation of a neutrino from r ~ to the outer edge of the Sun is adiabatic. That is, 
N e (r) changes sufficiently slowly that we may solve Schrodinger's equation for one r at 
a time, and then patch together the solutions. This means that our neutrino propagates 
outward through the Sun as one of the r-dependent eigenstates of the r-dependent Ti(r). 
Since the eigenvalues of 7i(r) do not cross at any r, and our neutrino is born in the 
heavier of the two r = eigenstates, it emerges from the Sun in the heavier of the two 
Tiy eigenstates [12]. The latter is the mass eigenstate we have called v% given according 
to Eq. (1.18) by 

^2 = ^e sin 6*0 + v x cos 0© . (1.23) 

Since this is an eigenstate of the vacuum Hamiltonian, the neutrino remains in it all the 
way to the surface of the Earth. The probability of observing the neutrino as a v e on 
Earth is then just the probability that v^ is a v e . That is [cf. Eq. (1.23)] [13], 

P (u e ^ u e ) = sin 2 9 Q . (1.24) 

We note that for 0© < 7r/4, this v e survival probability is less than 1/2. In contrast, when 
matter effects are negligible, the energy-averaged survival probability in two-neutrino 
oscillation cannot be less than 1/2 for any mixing angle [see Eq. (1.20)] [14]. 

II. The evidence for flavor metamorphosis, and what it has taught us: The 

persuasiveness of the evidence that neutrinos actually do change flavor in nature is 
summarized in Table 1.1. We discuss the different pieces of evidence, and what, together, 
they imply. 
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1. Neutrino mixing 7 



Table 1.1: The persuasiveness of the evidence for neutrino flavor change. The 
symbol L denotes the distance travelled by the neutrinos. LSND is the Liquid 
Scintillator Neutrino Detector experiment, and MiniBooNE is an experiment 
designed to confirm or refute LSND. 

Neutrinos Evidence for Flavor Change 

Atmospheric Compelling 

Accelerator (L = 250 and 735 km) Compelling 

Solar Compelling 

Reactor (L ~ 180 km) Compelling 

From Stopped fi + Decay (LSND) Unconfirmed by MiniBooNE 



The atmospheric neutrinos are produced in the Earth's atmosphere by cosmic rays, 
and then detected in an underground detector. The flux of cosmic rays that lead to 
neutrinos with energies above a few GeV is isotropic, so that these neutrinos are produced 
at the same rate all around the Earth. This can easily be shown to imply that at any 
underground site, the downward- and upward-going fluxes of multi-GeV neutrinos of a 
given flavor must be equal. That is, unless some mechanism changes the flux of neutrinos 
of the given flavor as they propagate, the flux coming down from zenith angle 9z must 
equal that coming up from angle tv — 9z [15]. 

The underground Super-Kamiokande (SK) detector finds that for multi-GeV 
atmospheric muon neutrinos, the 9z event distribution looks nothing like the expected 
6z <$ 7T — @z symmetric distribution. For cos 9z ~ 0.3, the observed v^ flux coming up 
from zenith angle tx — 9z is only about half that coming down from angle 9z [16]. Thus, 
some mechanism does change the v^ flux as the neutrinos travel to the detector. Since 
the upward-going muon neutrinos come from the atmosphere on the opposite side of the 
Earth from the detector, they travel much farther than the downward-going ones to reach 
the detector. Thus, if the muon neutrinos are oscillating away into another flavor, the 
upward-going ones have more distance (hence more time) in which to do so, which would 
explain why Flux Up < Flux Down. 

If atmospheric muon neutrinos are disappearing via oscillation into another flavor, then 
a significant fraction of accelerator-generated muon neutrinos should disappear on their 
way to a sufficiently distant detector. This disappearance has been observed by both the 
K2K [17] and MINOS [18] experiments. Each of these experiments measures its Vn flux in 
a detector near the neutrino source, before any oscillation is expected, and then measures 
it again in a detector 250 km from the source in the case of K2K, and 735 km from it in 
the case of MINOS. In its far detector, MINOS has observed 215 v^ events in a data 
sample where 336 ± 14.4 events would have been expected, in the absence of oscillation, 
on the basis of the near-detector measurements. Both K2K and MINOS also find that 
the energy spectrum of surviving muon neutrinos in the far detector is distorted in a way 
that is consistent with two-neutrino oscillation. 
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8 1 . Neutrino mixing 

The null results of short-baseline reactor neutrino experiments [8] imply limits on 
P(V^ — > PTj), which, assuming CPT invariance, are also limits on P(un — > v e ). From 
the latter, we know that the neutrinos into which the atmospheric, K2K, and MINOS 
muon neutrinos oscillate are not electron neutrinos, except possibly a small fraction of 
the time. All of the voluminous SK atmospheric neutrino data, corroborating data from 
other atmospheric neutrino experiments [19,20], K2K accelerator neutrino data, and 
existing MINOS accelerator neutrino data, are very well described by pure v^ — > v T 
quasi-two-neutrino oscillation. The allowed region for the oscillation parameters, Am| m 
and sin 2# a t m , which may be identified respectively with the parameters AM and 
4T /X (1 — Tn) in Eq. (1.15), is shown in Fig. 1.1. We note that this figure implies that at 
least one mass eigenstate Vi must have a mass exceeding 40 meV. 
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Figure 1.1: The region of the atmospheric oscillation parameters Am^ tm and 
sin 2 26» atm allowed by the SK, K2K, and MINOS data. The results of two different 
analyses of the SK ("Super K") data are shown [21]. 
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1. Neutrino mixing 9 

The neutrinos created in the Sun have been detected on Earth by several experiments, 
as discussed by K. Nakamura in this Review. The nuclear processes that power the Sun 
make only u e , not v^ or v T . For years, solar neutrino experiments had been finding that 
the solar v e flux arriving at the Earth is below the one expected from neutrino production 
calculations. Now, thanks especially to the Sudbury Neutrino Observatory (SNO), we 
have compelling evidence that the missing v e have simply changed into neutrinos of other 
flavors. 

SNO has studied the flux of high-energy solar neutrinos from B decay. This 
experiment detects these neutrinos via the reactions 

u + d^e~ +p + p , (1.25) 

v + d->v + p + n , (1.26) 

and 

v + e~ -»• u + e~ . (1.27) 

The first of these reactions, charged-current deuteron breakup, can be initiated only 
by a u e . Thus, it measures the flux 4>(u e ) of v e from B decay in the Sun. The 
second reaction, neutral-current deuteron breakup, can be initiated with equal cross 
sections by neutrinos of all active flavors. Thus, it measures <j>(v e ) + 4>( u h,t), where 
(/>(z/^ T ) is the flux of v^ and/or v T from the Sun. Finally, the third reaction, neutrino 
electron elastic scattering, can be triggered by a neutrino of any active flavor, but 
c(^,r e — »■ v^t e) — o~(v e e — ► v e e)/6.5. Thus, this reaction measures <f>(v e ) + 0(^ )T )/6.5. 
SNO finds from its observed rates for the two deuteron breakup reactions that [22] 



4>W, 



4> (Ve) + <P (vh,t) 



0.340 ± 0.023 (stat) +J;Jg (syst) . (1.28) 



Clearly, 0(^, r ) is not zero. This non-vanishing z/^ T flux from the Sun is "smoking-gun" 
evidence that some of the v e produced in the solar core do indeed change flavor. 

Corroborating information comes from the detection reaction ve~ — >■ ve~ , studied by 
both SNO and SK [23]. 

Change of neutrino flavor, whether in matter or vacuum, does not change the total 
neutrino flux. Thus, unless some of the solar v e are changing into sterile neutrinos, 
the total active high-energy flux measured by the neutral-current reaction (1.26) 
should agree with the predicted total 8 B solar neutrino flux based on calculations of 
neutrino production in the Sun. This predicted total is (5.49_ ' 81 ) x 10 6 cm _2 s _1 or 

(4.34_ ' 61 ) x 10 6 cm _2 s _1 , depending on assumptions about the solar heavy element 
abundances [24]. By comparison, the total active flux measured by reaction (1.26) 
is [4.94 ± 0.21 (stat) 1 ' 34 (syst)] x 10 6 cm _2 s _1 , in good agreement. This agreement 
provides evidence that neutrino production in the Sun is correctly understood, further 
strengthens the evidence that neutrinos really do change flavor, and strengthens the 
evidence that the previously-reported deficits of solar v e flux are due to this change of 
flavor. 
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10 1 . Neutrino mixing 

The strongly favored explanation of B solar neutrino flavor change is the LMA-MSW 
effect. As pointed out after Eq. (1.24), a v e survival probability below 1/2, which is 
indicated by Eq. (1.28), requires that solar matter effects play a significant role [25]. 
However, from Eq. (1.22) we see that as the energy E of a solar neutrino decreases, the 
vacuum (1st) term in the Hamiltonian 7i dominates more and more over the matter 
term. When we go from the B neutrinos with typical energies of ~6-7MeV to the 
monoenergetic Be neutrinos with energy 0.862 MeV, the matter term becomes fairly 
insignificant, and the v e survival probability is expected to be given by the vacuum 
oscillation formula of Eq. (1.20). In this formula, 9 is to be taken as the vacuum solar 
mixing angle 9q ~ 35° implied by the B solar neutrino data via Eqs. (1.28) and (1.24). 
When averaged over the energy-line shape, the oscillatory factor sin [1.27 Am (L/E)] is 
1/2, so that from Eq. (1.20) we expect that for the Be neutrinos, P [y e — > u e ) « 0.6. 

The Borexino experiment has now provided the first real time detection of the 
0.862 MeV Be solar neutrinos [26]. Borexino uses a liquid scintillator detector that 
detects these neutrinos via elastic neutrino-electron scattering. The experiment reports a 
'Be v e counting rate of [47 ± 7 (stat) ± 12(syst)] counts/day/lOOtons. Without any flavor 
change, this rate would have been expected to be [75 ±4] counts/day/lOOtons. With the 
degree of flavor change predicted by our understanding of the °B data [27] (see rough 
argument above), the rate would have been expected to be [49 ± 4] counts/day/lOOtons. 
The Borexino data are in nice agreement with the latter expectation, and the Borexino 
Collaboration is vigorously engaged in reducing its uncertainties. 

The LMA-MSW interpretation of B solar neutrino behavior implies that a substantial 
fraction of reactor V e that travel more than a hundred kilometers should disappear into 
antineutrinos of other flavors. The KamLAND experiment [28], which studies reactor 
V e that typically travel ~ 180 km to reach the detector, confirms this disappearance. 
In addition, KamLAND finds that the spectrum of the surviving V e that do reach the 
detector is distorted, relative to the no-oscillation spectrum. As Fig. 1.2 shows, the 
survival probability P (X^T — > X^r) measured by KamLAND is very well described by 
the hypothesis of neutrino oscillation. In particular, the measured survival probability 
displays the signature oscillatory behavior of the two-neutrino expression of Eq. (1.20). 
Ideally, the data in Fig. 1.2 would be plotted vs. L/E. However, KamLAND detects the 
V^ from a number of power reactors, at a variety of distances from the detector, so the 
distance L travelled by any given X^T is unknown. Consequently, Fig. 1.2 plots the data 
vs. Lq/E, where Lq = 180 km is a flux- weighted average travel distance. The oscillation 
curve and histogram in the figure take the actual distances to the individual reactors 
into account. Nevertheless, almost two cycles of the sinusoidal structure expected from 
neutrino oscillation are still plainly visible. 

The region allowed by solar neutrino experiments for the two-neutrino vacuum 
oscillation parameters Am Q and #©, and that allowed by KamLAND for what we believe 
to be the same parameters, are shown in Fig. 1.3. From this figure, we see that there is a 
region of overlap. This is strong evidence that the behavior of both solar neutrinos and 
reactor antineutrinos has been correctly understood. A joint analysis of KamLAND and 
solar neutrino data assuming CPT invariance yields AttIq = (7.59 ± 0.21) x 10 eV 
and tan 2 0© = 0.47+{J;[jf [28]. 
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Figure 1.2: Ratio of the background- and geo-neutrino subtracted v e spectrum to 
the no-oscillation expectation as a function of Lq/E [28]. See text for explanation. 

That 6> a t m and 9@ are both large, in striking contrast to all quark mixing angles, is 
very interesting. 

The neutrinos studied by the LSND experiment [29] come from the decay fi + —> e~^u e Un 
of muons at rest. While this decay does not produce V e , an excess of V e over expected 
background is reported by the experiment. This excess is interpreted as due to oscillation 
of some of the Vn produced by fx + decay into V e . The related Karlsruhe Rutherford 
Medium Energy Neutrino (KARMEN) experiment [30] sees no indication for such an 
oscillation. However, the LSND and KARMEN experiments are not identical; at LSND 
the neutrino travels a distance L « 30 m before detection, while at KARMEN it travels 
L « 18 m. The KARMEN results exclude a portion of the neutrino parameter region 
favored by LSND, but not all of it. A joint analysis [31] of the results of both experiments 
finds that a splitting 0.2 < Am^g^ ^ leV 2 and mixing 0.003 ^ sin 2 26*lsnd ~ 0.03, 
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Figure 1.3: Regions allowed by the "solar" neutrino oscillation parameters by 
KamLAND and by solar neutrino experiments [28]. 



or a splitting Am^ SND 
experiments. 



7eV and mixing sin 2#l§nd ~ 0.004, might explain both 



To confirm or exclude the LSND oscillation signal, the MiniBooNE experiment was 
launched. MiniBooNE studies Vu and VZ that travel a distance L of 540 m and have a 
typical energy E of 700 MeV, so that L/E is of order 1 km/GeV as in LSND. MiniBooNE's 
first results [32], regarding a search for i/n — > v e oscillation in a Vn beam, do not confirm 
LSND. For neutrino energies 475 < E < 3000 MeV, there is no significant excess of events 
above background. A joint analysis of the MiniBooNE data at these energies and the 



LSND data excludes at 98% CL two-neutrino v7l—*v^ oscillation as an explanation of the 
LSND FjT excess. To be sure, there is an excess of MiniBooNE v e candidate events below 
475 MeV. This low-energy excess cannot be explained by two-neutrino oscillation, and 
its source is being studied. Possibilities include an unidentified background, a Standard 
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Model effect that has been proposed only recently [33], and many-neutrino oscillation 
with a CP violation that allows the antineutrino oscillation reported by LSND to differ 
from the neutrino results reported so far by MiniBooNE [34] . 

The MiniBooNE detector is illuminated by both the neutrino beam constructed for 
the purpose, and the beam that is aimed at the MINOS detector. The distance L to 
MiniBooNE from the neutrino source is 40% larger in the latter beam than in the former. 
When matter effects may be neglected, the probability of oscillation depends on L and 
the beam energy E only through L/E [cf. Eq. (1.9)]. Thus, if the low-energy excess seen 
by MiniBooNE is neutrino oscillation, it should appear at a 40% higher energy in the 
beam directed at MINOS than in MiniBooNE's own beam. Whether it does or not is 
under investigation. 

The regions of neutrino parameter space favored or excluded by various neutrino 
oscillation experiments are shown in Fig. 1.4. 

III. Neutrino spectra and mixings: If there are only three neutrino mass eigenstates, 
v\, h>2, and z/3, then there are only three mass splittings AmL and they obviously satisfy 

Am^ + Am|i + Amf 3 = . (1.29) 

However, as we have seen, the Am values required to explain the flavor changes of the 
atmospheric, solar, and LSND neutrinos are of three different orders of magnitude. Thus, 
they cannot possibly obey the constraint of Eq. (1.29). If all of the reported changes of 
flavor are genuine, then nature must contain at least four neutrino mass eigenstates [35]. 
As explained in Sec. I, one linear combination of these mass eigenstates would have to be 
sterile. 

If further MiniBooNE results do not confirm the LSND oscillation, then nature 
may well contain only three neutrino mass eigenstates. The neutrino spectrum then 
contains two mass eigenstates separated by the splitting Attiq needed to explain the 
solar and KamLAND data, and a third eigenstate separated from the first two by the 
larger splitting Am^ tm called for by the atmospheric, MINOS, and K2K data. Current 
experiments do not tell us whether the solar pair - - the two eigenstates separated 
by AttIq - - is at the bottom or the top of the spectrum. These two possibilities are 
usually referred to, respectively, as a normal and an inverted spectrum. The study of 
flavor changes of accelerator-generated neutrinos and antineutrinos that pass through 
matter can discriminate between these two spectra (see Sec. V). If the solar pair is at 
the bottom, then the spectrum is of the form shown in Fig. 1.5. There we include the 
approximate flavor content of each mass eigenstate, the flavor-cn fraction of eigenstate v$ 
being simply |(i/ a |i^)| = | C^cti I - The flavor content shown assumes that the atmospheric 
mixing angle is maximal, which gives the best fit to the atmospheric data [16] and, as 
indicated in Fig. 1.1, to the MINOS data. The content shown also takes into account 
the now-established LMA-MSW explanation of solar neutrino behavior. For simplicity, it 
neglects the small, as- yet- unknown z/ e fraction of z/3 (see below) . 

When there are only three neutrino mass eigenstates, and the corresponding three 
familiar neutrinos of definite flavor, the leptonic mixing matrix U can be written as 
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Figure 1.4: The regions of squared- mass splitting and mixing angle favored or 
excluded by various experiments. This figure was contributed by H. Murayama 
(University of California, Berkeley). References to the data used in the figure can 
be found at http://hitoshi.berkeley.edu/neutrino/. 
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Figure 1.5: A three-neutrino squared-mass spectrum that accounts for the observed 
flavor changes of solar, reactor, atmospheric, and long-baseline accelerator neutrinos. 
The v e fraction of each mass eigenstate is crosshatched, the Vn fraction is indicated 
by right-leaning hatching, and the v T fraction by left-leaning hatching. 



U 





v\ 


Ve 


C12C13 


Vfi 


-S12C23 - Ci2S23«13e 


V T 


^12^23 - Ci 2 C 2 3Sl3e 



x diag ( e 



,«*i/2 



id 
id 

j/2 



V2 

S12C13 

C12C23 - Sl2S23Sl3e 

-C12«23 - Sl2C23«13e 



iS 
iS 



V'i 

«l3e 
S23C13 
C23C13 



iS 



(1.30) 



Here, v\ and v^ are the members of the solar pair, with m<i > mi, and v% is the 
isolated neutrino, which may be heavier or lighter than the solar pair. Inside the matrix, 
Cij = cos 9ij and Sij = sinO™, where the three O^s are mixing angles. The quantities 
5, a>i, and «2 are CP-violating phases. The phases a\ and «2, known as Majorana 
phases, have physical consequences only if neutrinos are Majorana particles, identical 
to their antiparticles. Then these phases influence neutrinoless double-beta decay [see 
Sec. IV] and other processes [36]. However, as we see from Eq. (1.9), a\ and «2 do 
not affect neutrino oscillation, regardless of whether neutrinos are Majorana particles. 
Apart from the phases ai, «2, which have no quark analogues, the parametrization of 
the leptonic mixing matrix in Eq. (1.30) is identical to that [37] advocated for the quark 
mixing matrix by Ceccucci, Ligeti, and Sakai in their article in this Review. 

From bounds on the short-distance oscillation of reactor V e [8] and other data, at 
2 a, \U e ^\ 2 ^ 0.032 [38]. (Thus, the v e fraction of ^3 would have been too small to see in 
Fig. 1.5; this is the reason it was neglected.) From Eq. (1.30), we see that the bound on 
I C/e3 1 implies that s\^ ^ 0.032. From Eq. (1.30), we also see that the CP-violating phase 
5, which is the sole phase in the U matrix that can produce CP violation in neutrino 
oscillation, enters U only in combination with S13. Thus, the size of CP violation in 
oscillation will depend on S13. 
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Given that S13 is small, Eqs. (1.30), (1.15), and (1.17) imply that the atmospheric 
mixing angle a tm extracted from v^ disappearance measurements is approximately 6*23, 
while Eqs. (1.30) and (1.18) (with v a = v e and 9 = 6@) imply that 6>© ~ 6*12- 

IV. The neutrino- antineutrino relation: Unlike quarks and charged leptons, 
neutrinos may be their own antiparticles. Whether they are depends on the nature of the 
physics that gives them mass. 

In the Standard Model (SM), neutrinos are assumed to be massless. Now that we 
know they do have masses, it is straightforward to extend the SM to accommodate these 
masses in the same way that this model accommodates quark and charged lepton masses. 
When a neutrino v is assumed to be massless, the SM does not contain the chirally 
right-handed neutrino field v R , but only the left-handed field v^ that couples to the W 
and Z bosons. To accommodate the v mass in the same manner as quark masses are 
accommodated, we add v R to the Model. Then we may construct the "Dirac mass term" 

C D = -m D V L v R + h.c. , (1.31) 

in which mp is a constant. This term, which mimics the mass terms of quarks and 
charged leptons, conserves the lepton number L that distinguishes neutrinos and 
negatively- charged leptons on the one hand from antineutrinos and positively-charged 
leptons on the other. Since everything else in the SM also conserves L, we then have an 
L-conserving world. In such a world, each neutrino mass eigenstate z/j differs from its 
antiparticle V^, the difference being that L(Vi) = —L{yj). When V{ 7^ vi, we refer to the 
V{ — V{ complex as a "Dirac neutrino." 

Once v R has been added to our description of neutrinos, a "Majorana mass term," 

C M = -m R v c R v R + h.c. , (1.32) 

can be constructed out of v R and its charge conjugate, v R . In this term, m R is another 
constant. Since both v R and v R absorb v and create V, Cm mixes v and V. Thus, a 
Majorana mass term doe s not conserve L. In somewhat the same way that, neglecting 
CP violation, K — K Q m ixin g causes the neutral kaon mass eigenstates to be the 
self-conjugate states (K ± iir )/v / 2, the v — v mixing induced by a Majorana mass term 
causes the neutrino mass eigenstates to be self-conjugate: T>1 — v^. That is, for a given 
helicity h, Vi{h) = vi{h). We then refer to V{ as a "Majorana neutrino." 

Suppose the right-handed neutrinos required by Dirac mass terms have been added to 
the SM. If we insist that this extended SM conserve L, then, of course, Majorana mass 
terms are forbidden. However, if we do not impose L conservation, but require only the 
general principles of gauge invariance and renormalizability, then Majorana mass terms 
like that of Eq. (1.32) are expected to be present. As a result, L is violated, and neutrinos 
are Majorana particles [39]. 

In the see-saw mechanism [40] , which is the most popular explanation of why neutrinos 
- although massive -- are nevertheless so light, both Dirac and Majorana mass terms 
are present. Hence, the neutrinos are Majorana particles. However, while half of them 
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are the familiar light neutrinos, the other half are extremely heavy Majorana particles 
referred to as the Ni, with masses possibly as large as the GUT scale. The iVj may have 
played a crucial role in baryogenesis in the early universe, as we shall discuss in Sec. V. 

How can the theoretical expectation that nature contains Majorana mass terms, so 
that L is violated and neutrinos are Majorana particles, be confirmed experimentally? 
The promising approach is to search for neutrinoless double-beta decay (0vfif3). This 
is the process (A, Z) — > (A, Z + 2) + 2e~, in which a nucleus containing A nucleons, 
Z of which are protons, decays to a nucleus containing Z + 2 protons by emitting two 
electrons. While 0^/3/3 can in principle receive contributions from a variety of mechanisms 
(R-parity- violating supersymmetric couplings, for example), it is easy to show explicitly 
that its observation at any non-vanishing rate would imply that nature contains at least 
one Majorana neutrino mass term [41]. The neutrino mass eigenstates must then be 
Majorana neutrinos. 

Quarks and charged leptons cannot have Majorana mass terms, because such terms 
mix fermion and antifermion, and q «-> q or £ ^ £ would not conserve electric charge. 
Thus, the discovery of Qv{3(3 would demonstrate that the physics of neutrino masses is 
unlike that of the masses of all other fermions. 

The dominant mechanism for Ozv/3/3 is expected to be the one depicted in Fig. 1.6. 
There, a pair of virtual W bosons are emitted by the parent nucleus, and then these 
W bosons exchange one or another of the light neutrino mass eigenstates Vi to produce 
the outgoing electrons. The 0^/3/3 amplitude is then a sum over the contributions of the 
different v^. It is assumed that the interactions at the two leptonic W vertices are those 
of the SM. 
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Figure 1.6: The dominant mechanism for 0v(3/3. The diagram does not exist unless 



Va = Vi. 



Since the exchanged v^ is created together with an e~ , the left-handed SM current that 
creates it gives it the helicity we associate, in common parlance, with an "antineutrino." 
That is, the V{ is almost totally right-handed, but has a small left-handed-helicity 
component, whose amplitude is of order rrii/E, where E is the v^ energy. At the vertex 
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where this V{ is absorbed, the absorbing left-handed SM current can absorb only its 
small left-handed-helicity component without further suppression. Consequently, the 
z/j-exchange contribution to the 0^/3/3 amplitude is proportional to m^. From Fig. 1.6, 
we see that this contribution is also proportional to L/ 2 -. Thus, summing over the 
contributions of all the z/j, we conclude that the amplitude for 0vf3f3 is proportional to the 
quantity 



Yl mi U ei 



< mpp > | , (1.33) 



commonly referred to as the "effective Majorana mass for neutrinoless double-beta 
decay" [42]. 

To how small an | < man > | should a 0z//3/9 search be sensitive? In answering this 
question, it makes sense to assume there are only three neutrino mass eigenstates - - if 
there are more, | < m^ > | might be larger. Suppose that there are just three mass 
eigenstates, and that the solar pair, ui and z/2, is at the top of the spectrum, so that we 
have an inverted spectrum. If the various v^ are not much heavier than demanded by the 
observed splittings Am 2 tm and Am , then in | < mag > |, Eq. (1.33), the contribution 
of z/3 may be neglected, because both 7713 and \U^\ = s 2 3 are small. From Eqs. (1.33) 
and (1.30), approximating C13 by unity, we then have that 



I < mpp> I ~m o Wl-sin 2 26> sin 2 f-^j . (1.34) 

Here, mo is the average mass of the members of the solar pair, whose splitting will be 
invisible in a practical 0u/3/3 experiment, and Act = «2 — ol\ is a CP-violating phase. 
Although Act is completely unknown, we see from Eq. (1.34) that 

I < m (3(3 > I — m,QCOs29Q . (1.35) 



Now, in an inverted spectrum, tjiq > wAm 2 tm . At 90% CL, J Am 2 tm > 45meV [see 

Fig. 1.1], while at 95% CL, cos 26*0 > 0.25 [see Fig. 1.3]. Thus, if neutrinos are Majorana 
particles, and the spectrum is as we have assumed, a Ovflfl experiment sensitive to 
I < m f3/3 > I ~ 10 meV would have an excellent chance of observing a signal. If the 
spectrum is inverted, but the vi masses are larger than the Am 2 tm - and AmQ-demanded 
minimum values we have assumed above, then once again | < mpp > \ is larger than 10 
meV [43], and an experiment sensitive to 10 meV still has an excellent chance of seeing a 
signal. 

If the solar pair is at the bottom of the spectrum, rather than at the top, then 
I < m /3/3 > I is n °t as tightly constrained, and can be anywhere from the present bound 
of 0.3-1.0 eV down to invisibly small [43,44]. For a discussion of the present bounds, see 
the article by Vogel and Piepke in this Review [45]. 

V. Questions to be answered: The strong evidence for neutrino flavor metamorphosis 
- hence neutrino mass -- opens many questions about the neutrinos. These questions, 
which hopefully will be answered by future experiments, include the following: 
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i) How many neutrino species are there? Do sterile neutrinos exist? 

This question is being addressed by the MiniBooNE experiment [32]. If MiniBooNE's 
final result is positive, the implications will be far-reaching. We will have learned that 
either there are more than three neutrino species and at least one of these species is 
sterile, or else there is an even more amazing departure from what has been our picture 
of the neutrino world. 

ii) What are the masses of the mass eigenstates vi ? 

Assuming there are only three z^, we need to find out whether the solar pair, z/12, is at 
the bottom of the spectrum or at its top. This can be done by exploiting matter effects in 
long-baseline neutrino and antineutrino oscillations. These matter effects will determine 
the sign one wishes to learn -- that of {777.3 — t( m 2 ^~ m i)/2]} ~~ relative to a sign that 
is already known -- that of the interaction energy of Eq. (1.21). Grand unified theories 
favor a spectrum with the closely spaced solar pair at the bottom [46]. The neutrino 
spectrum would then resemble the spectra of the quarks, to which grand unified theories 
relate the neutrinos. A neutrino spectrum with the closely spaced solar pair at the top 
would be quite un-quark-like, and would suggest the existence of a new symmetry that 
leads to the near degeneracy at the top of the spectrum. 

While flavor-change experiments can determine a spectral pattern such as the one 
in Fig. 1.5, they cannot tell us the distance of the entire pattern from the zero of 
squared-mass. One might discover that distance via study of the j3 energy spectrum 
in tritium (3 decay, if the mass of some V{ with appreciable coupling to an electron is 
large enough to be within reach of a feasible experiment. One might also gain some 
information on the distance from zero by measuring | < 777/3/3 > 1 > the effective Majorana 
mass for neutrinoless double-beta decay [43-45] (see Vogel and Piepke in this Review). 
Finally, one might obtain information on this distance from cosmology or astrophysics. 
Indeed, from current cosmological data and some cosmological assumptions, it is already 
concluded that [47] 

"mj < (0.17-2.0) eV . (1.36) 



E 



Here, the sum runs over the masses of all the light neutrino mass eigenstates v^ that may 
exist and that were in thermal equilibrium in the early universe. The range quoted in 
Eq. (1.36) reflects the dependence of this upper bound on the underlying cosmological 
assumptions and on which data are used [47]. 

If there are just three t/j, and their spectrum is either the one shown in Fig. 1.5 or its 
inverted version, then Eq. (1.36) implies that the mass of the heaviest i/j, Mass [Heaviest 
z/j], cannot exceed (0.07 - 0.7) eV. Moreover, Mass [Heaviest i/j] obviously cannot be less 



than J Am^ tm , which in turn is not less than 0.04 eV, as previously noted. Thus, if the 
cosmological assumptions behind Eq. (1.36) are correct, then 

0.04 eV < Mass [Heaviest v t ] < (0.07 - 0.7) eV . (1.37) 
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Hi) Are the neutrino mass eigenstates Majorana particles ? 

The confirmed observation of neutrinoless double-beta decay would establish that the 
answer is "yes." If there are only three z/j, knowledge that the spectrum is inverted and 
a definitive upper bound on | < mpg > | that is well below 0.01 eV would establish 
(barring exotic contributions to Oz/,9/3) that the answer is "no" [see discussion after 

Eq. (1.35)] [43,44]. 

iv) What are the mixing angles in the leptonic mixing matrix U ? 

The solar mixing angle 9@ ~ #12 is already rather well determined. 

The atmospheric mixing angle # a tm — #23 is constrained by the most stringent analysis 
to lie, at 90% CL, in the region where sin 2# a t m > 0.92 [16]. This region is still fairly 
large: 37° to 53°. A more precise value of sin 26* a tm 5 and, in particular, its deviation 
from unity, can be sought in precision long-baseline v^ disappearance experiments. If 
sin 2# a tm 7^ 1, so that 6* a t m 7^ 45°, one can determine whether it lies below or above 45° 
with the help of a reactor V e experiment [48,49]. Once we know whether the neutrino 
spectrum is normal or inverted, this determination will tell us whether the heaviest mass 
eigenstate is more v T than z/^, as naively expected, or more v^ than v T [cf. Eq. (1.30)]. 

A knowledge of the small mixing angle #13 is important not only to help complete our 
picture of leptonic mixing, but also because, as Eq. (1.30) made clear, all CP-violating 
effects of the phase 5 are proportional to sin 6*13. Thus, a knowledge of the order of 
magnitude of #13 would help guide the planning of experiments to probe CP violation. 
From Eq. (1.30), we recall that sin #13 is the v e fraction of z/3. The z/3 is the isolated 
neutrino that lies at one end of the atmospheric squared-mass gap Am^ tm , so an 
experiment seeking to measure #13 should have an L/E that makes it sensitive to 
Am| m , and should involve z/ e . Planned approaches include a sensitive search for the 
disappearance of reactor z/ e while they travel a distance L ~ 1 km, and an accelerator 
neutrino search for z/^ — > z/ e and VZ, — > V^ with a beamline L > several hundred km. 

If LSND is confirmed, then (barring the still more revolutionary) the matrix U is at 
least 4x4, and contains many more than three angles. A rich program, including short 
baseline experiments with multiple detectors, will be needed to learn about both the 
squared-mass spectrum and the mixing matrix. 

Given the large sizes of # a tm and #0, we already know that leptonic mixing is 
very different from its quark counterpart, where all the mixing angles are small. This 
difference, and the striking contrast between the tiny neutrino masses and the very much 
larger quark masses, suggest that the physics underlying neutrino masses and mixing may 
be very different from the physics behind quark masses and mixing. 

v) Does the behavior of neutrinos violate CP? 

From Eqs. (1.9), (1.13), and (1.30), we see that if the CP-violating phase 5 
and the small mixing angle #13 are both non-vanishing, there will be CP-violating 
differences between neutrino and antineutrino oscillation probabilities. Observation of 
these differences would establish that CP violation is not a peculiarity of quarks. 

The CP-violating difference P(v a — ► up) — P{v a — > Vp) between "neutrino" and 
"antineutrino" oscillation probabilities is independent of whether the mass eigenstates 
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V{ are Majorana or Dirac particles. To study v^ — > v e with a super-intense but 
conventionally generated neutrino beam, for example, one would create the beam via 
the process tv + — ► /i + v^ and detect it via v^ + target — > e~ + . . .. To study V^ — ► z7 e , 
one would create the beam via 7r~ — ► |U _ 77^, and detect it via V^ + target — ► e + + . . .. 
Whether 77^ = v^ or not, the amplitudes for the latter two processes are proportional to 
U^i and U*i, respectively. In contrast, the amplitudes for their v^ — > v e counterparts 
are proportional to U*^ and XJ e {. As this illustrates, Eq. (1.13) relates "neutrino" and 
"antineutrino" oscillation probabilities even when the neutrino mass eigenstates are their 
own antiparticles. 

The baryon asymmetry of the universe could not have developed without some violation 
of CP during the universe's early history. The one known source of CP violation — the 
complex phase in the quark mixing matrix -- could not have produced sufficiently large 
effects. Thus, perhaps leptonic CP violation is responsible for the baryon asymmetry. 
The see-saw mechanism predicts very heavy Majorana neutral leptons Nj, (see Sec. IV), 
which would have been produced in the Big Bang. Perhaps CP violation in the leptonic 
decays of an N$ led to the inequality 

r (Ni - £+ + ...) ± r (Ni - r + .. .) , (1.38) 

which would have resulted in unequal numbers of £ + and t~ in the early universe [50]. 
This leptogenesis could have been followed by nonperturbative SM processes that would 
have converted the lepton asymmetry, in part, into the observed baryon asymmetry [51]. 

While the connection between the CP violation that would have led to leptogenesis, 
and that which we hope to observe in neutrino oscillation, is model-dependent, it is 
not likely that we have either of these without the other [52], because in the see-saw 
picture, these two CP violations both arise from the same matrix of coupling constants. 
This makes the search for CP violation in neutrino oscillation very interesting indeed. 
Depending on the rough size of #13, this CP violation may be observable with a very 
intense conventional neutrino beam, or may require a "neutrino factory," whose neutrinos 
come from the decay of stored muons or radioactive nuclei. The detailed study of CP 
violation may require a neutrino factory in any case. 

With a conventional beam, one would seek CP violation, and try to determine whether 
the mass spectrum is normal or inverted, by studying the oscillations v^ — > v e and 
z7^ — * Pe- The appearance probability for v e in a beam that is initially Vu can be written 
for sin 2 2#i 3 < 0.2 [53] 

P{^ -> v e ) = sin 2 2#i3 Ti -a sin 2#i 3 T 2 + a sin 29 13 T 3 + a 2 T A . (1.39) 

Here, a = Am^/Am^ is the small (~ 1/30) ratio between the solar and atmospheric 
squared-mass splittings, and 

9 sin 2 [(1 -x)A] . 

Ti=sin 2 2 3 77 T^- 1 , (1-40) 

(1-x) 

r^ <■ a sin (xA) sin [(1 — x) A] , , . 

T 2 = sin 5 sin 2^12 sin 2^3 sin A i '- [ ). {—!■ , 1.41 

x (1 — X) 
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r-r, <- . ~„ . ~„ a sin (xA) sin [(1 — x) A] ._ ,_ 

T 3 = cos 5 sin 2^2 sin 2^3 cos A i ^ £ -^ , (1.42) 

x (1 — X) 

and 

r 4 = cos 2 %sin 2 2^ 12 Sin2( 9 XA) . (1.43) 

x z 

In these expressions, A = Am^L/AE is the kinematical phase of the oscillation. The 
quantity x = 2\/2GpN e E / Am^, with Gp the Fermi coupling constant and N e the 
electron number density, is a measure of the importance of the matter effect resulting 
from coherent forward-scattering of electron neutrinos from ambient electrons as the 
neutrinos travel through the earth from the source to the detector [cf. Sec. I]. In 
the appearance probability P(^ —* v e ), the T\ term represents the oscillation due 
to the atmospheric-mass-splitting scale, the T4 term represents the oscillation due 
to the solar-mass-splitting scale, and the T2 and T3 terms are the CP-violating and 
CP-conserving interference terms, respectively. 

The probability for the corresponding antineutrino oscillation, P(l^7 — *• T^), is the same 
as the probability P(^ — *■ v e ) given by Eqs. (1.39)-(1.43), but with the signs in front of 
both x and sin 6 reversed: both the matter effect and CP violation lead to a difference 
between the v^ — > u e and 1^7 — ^ ^e oscillation probabilities. In view of the dependence of 
x on Am^ , and in particular on the sign of Am^ , the matter effect can reveal whether 
the neutrino mass spectrum is normal or inverted. However, to determine the nature 
of the spectrum, and to establish the presence of CP violation, it obviously will be 
necessary to disentangle the matter effect from CP violation in the neutrino-antineutrino 
oscillation probability difference that is actually observed. To this end, complementary 
measurements will be extremely important. These can take advantage of the differing 
dependences on the matter effect and on CP violation in P(z^ — > u e ). 

vi) Will we encounter the completely unexpected? 

The study of neutrinos has been characterized by surprises. It would be surprising if 
further surprises were not in store. The possibilities include new, non-Standard-Model 
interactions, unexpectedly large magnetic and electric dipole moments [54] , unexpectedly 
short lifetimes, and violations of CPT invariance, Lorentz invariance, or the equivalence 
principle. 

The questions we have discussed, and other questions about the world of neutrinos, 
will be the focus of a major experimental program in the years to come. 

Acknowledgements 

I am grateful to Susan Kayser for her crucial role in the production of this manuscript. 

References: 

1. This matrix is sometimes referred to as the Maki-Nakagawa- Sakata matrix, or 
as the Pontecorvo-Maki-Nakagawa-Sakata matrix, in recognition of the pioneering 
contributions of these scientists to the physics of mixing and oscillation. See Z. Maki, 
M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962); 

B. Pontecorvo, Zh. Eksp. Teor. Fiz. 53, 1717 (1967) [Sov. Phys. JETP 26, 984 

(1968)]. 

July 14, 2008 16:26 



1. Neutrino mixing 23 

2. D. Karlen in this Review. 

3. B. Kayser, Phys. Rev. D24, 110 (1981); 

F. Boehm and P. Vogel, Physics of Massive Neutrinos (Cambridge University Press, 
Cambridge, 1987) p. 87; 

C. Giunti, C. Kim, and U. Lee, Phys. Rev. D44, 3635 (1991); 

J. Rich, Phys. Rev. D48, 4318 (1993); 

H. Lipkin, Phys. Lett. B348, 604 (1995); 

W. Grimus and P. Stockinger, Phys. Rev. D54, 3414 (1996); 

T. Goldman, hep-ph/9604357; 

Y. Grossman and H. Lipkin, Phys. Rev. D55, 2760 (1997); 

W. Grimus, S. Mohanty, and P. Stockinger, in Proc. of the 17th Int. Workshop 
on Weak Interactions and Neutrinos, eds. C. Dominguez and R. Viollier (World 
Scientific, Singapore, 2000) p. 355; 

L. Stodolsky, Phys. Rev. D58, 036006 (1998); 

C. Giunti, Phys. Scripta 67, 29 (2003); 

M. Beuthe, Phys. Rept. 375, 105 (2003) and Phys. Rev. D66, 013003 (2002), and 
references therein; 

H. Lipkin, Phys. Lett. B579, 355 (2004); 

H. Lipkin, Phys. Lett. B642, 366 (2006); 

C. Giunti and C. Kim, Fundamentals of Neutrino Physics and Astrophysics (Oxford 
University Press, Oxford, 2007), p. 283. 

4. H. Lipkin (2006), Ref. 3. 

5. G. Fogli, E. Lisi, and G. Scioscia, Phys. Rev. D52, 5334 (1995). 

6. L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

7. S. Mikheyev and A. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 913 
(1986)]; Zh. Eksp. Teor. Fiz. 91, 7, (1986) [Sov. Phys. JETP 64, 4 (1986)]; Nuovo 
Cimento 9C, 17 (1986). 

8. The Bugey Collaboration (B. Achkar et al), Nucl. Phys. B434, 503 (1995); 

The Palo Verde Collaboration (F. Boehm et al), Phys. Rev. D64, 112001 (2001); 
The CHOOZ Collaboration (M. Apollonio et al), Eur. Phys. J. C27, 331 (2003). 

9. P. Langacker, J. Leveille, and J. Sheiman, Phys. Rev. D27, 1228 (1983); 
The corresponding energy for anti-neutrinos is negative. 

10. G. L. Fogli, E. Lisi, and D. Montanino, Phys. Rev. D54, 2048 (1996); 

A. de Gouvea, A. Friedland, and H. Murayama, Phys. Lett. B490, 125 (2000). 

11. J. Bahcall, Neutrino Astrophysics, (Cambridge Univ. Press, Cambridge, UK 1989). 

12. A more quantitative analysis has shown that this is true to better than 90% 
probability. See H. Nunokawa, S. Parke, and R. Zukanovich Funchal, Phys. Rev. 
D74, 013006 (2006). 

July 14, 2008 16:26 



24 1 . Neutrino mixing 

13. S. Parke, Phys. Rev. Lett. 57, 1275 (1986). 

14. We thank J. Beacom and A. Smirnov for invaluable conversations on how LMA-MSW 
works. For an early description, see S. Mikheyev and A. Smirnov, Ref. 7 (first paper). 

15. D. Ayres et al, in Proc. of the 1982 DPF Summer Study on Elementary Particle 
Physics and Future Facilities, p. 590; 

G. Dass and K. Sarma, Phys. Rev. D30, 80 (1984); 

J. Flanagan, J. Learned, and S. Pakvasa, Phys. Rev. D57, 2649 (1998); 

B. Kayser, in Proc. of the 17th Int. Workshop on Weak Interactions and Neutrinos, 
eds. C. Dominguez and R. Viollier (World Scientific, Singapore, 2000) p. 339. 

16. The Super-Kamiokande Collaboration (Y. Ashie et a/.), Phys. Rev. D71, 112005 
(2005). 

17. The K2K Collaboration (M. Ahn et a/.), Phys. Rev. D74, 072003 (2006). 

18. The MINOS Collaboration (D. Michael et al), Phys. Rev. Lett. 97, 191801 (2006). 

19. MACRO Collaboration (M. Ambrosio et al), Phys. Lett. B566, 35 (2003); 
MACRO Collaboration, Eur. Phys. J. C, 36, 323 (2004). 

20. M. Sanchez et al., Phys. Rev. D68, 113004 (2003); 
W.W.M. Allison et al., Phys. Rev. D72, 052005 (2005). 

21. This figure is taken from a talk given by Niki Saoulidou on behalf of the MINOS 
Collaboration on July 19, 2007. The Super-Kamiokande contours shown are based on 
Ref. 16 and on The Super-Kamiokande Collaboration (Y. Ashie et al.), Phys. Rev. 
Lett. 93, 101801 (2004). 

22. The SNO Collaboration (B. Aharmim et al.), Phys. Rev. C72, 055502 (2005). 

23. Y. Koshio, in 2003 Electroweak Interactions and Unified Theories (Proceedings of the 
38th Rencontres de Moriond), ed. J. Tran Thanh Van (The Gioi, Vietnam, 2003) 
p. 3; 

The Super-Kamiokande Collaboration (J. Hosaka et al.), Phys. Rev. D73, 112001 
(2006). 

24. J. Bahcall, S. Basu, and A. Serenelli, Astrophys. J. Suppl. 165, 400 (2006). 

25. G. Fogli et al, Phys. Lett. B583, 149 (2004). 

26. The Borexino Collaboration (C. Arpesella et al), Phys. Lett. B658, 101 (2008). 

27. J. Bahcall, M.C. Gonzalez-Garcia, and C. Pena-Garay, JHEP 0408, 16 (2004); 
G. Fogli et al, Prog. Part. Nucl. Phys. 57, 742 (2006). 

28. The KamLAND Collaboration (S. Abe et al), arXiv: 0801.4589. 

29. The LSND Collaboration (A. Aguilar et al), Phys. Rev. D64, 112007 (2001). 

30. The KARMEN Collaboration (B. Armbruster et al), Phys. Rev. D65, 112001 
(2002). 

31. E. Church et al, Phys. Rev. D66, 013001 (2003). 

32. The MiniBooNE Collaboration (A. Aguilar- Arevalo et al), Phys. Rev. Lett. 98, 
231801 (2007). 

July 14, 2008 16:26 



1. Neutrino mixing 25 

33. J. Harvey, C. Hill, and R. Hill, Phys. Rev. Lett. 99, 261601 (2007). 

34. G. Karagiorgi et al, Phys. Rev. D75, 013011 (2007); 

See, however, M. Maltoni and T. Schwetz, Phys. Rev. D76, 093005 (2007). 

35. For an alternative possibility entailing CPT violation, see H. Murayama and T. 
Yanagida, Phys. Lett. B520, 263 (2001); 

G. Barenboim et al, JHEP 0210, 001 (2002); 

However, after KamLAND, this alternative is disfavored. M. C. Gonzalez-Garcia, M. 
Maltoni, and T. Schwetz, Phys. Rev. D68, 053007 (2003); 

G. Barenboim, L. Borissov, and J. Lykken, hep-ph/0212116 v2. 

36. J. Schechter and J. Valle, Phys. Rev. D23, 1666 (1981); 
J. Nieves and P. Pal, Phys. Rev. D64, 076005 (2001); 

A. de Gouvea, B. Kayser, and R. Mohapatra, Phys. Rev. D67, 053004 (2003). 

37. L.-L. Chau and W.-Y. Keung, Phys. Rev. Lett. 53, 1802 (1984); 
H. Harari and M. Leurer, Phys. Lett. B181, 123 (1986); 

F.J. Botella and L.-L. Chau, Phys. Lett. B168, 97 (1986); 
H. Fritzsch and J. Plankl, Phys. Rev. D35, 1732 (1987). 

38. G. Fogli et al, see Ref. 27. 

39. We thank Belen Gavela for introducing us to this argument. 

40. M. Gell-Mann, P. Ramond, and R. Slansky, in: Supergravity, eds. D. Freedman and 
P. van Nieuwenhuizen (North Holland, Amsterdam, 1979) p. 315; 

T. Yanagida, in: Proceedings of the Workshop on Unified Theory and Baryon Number 
in the Universe, eds. O. Sawada and A. Sugamoto (KEK, Tsukuba, Japan, 1979); 

R. Mohapatra and G. Senjanovic: Phys. Rev. Lett. 44, 912 (1980) and Phys. Rev. 
D23, 165 (1981); 

P. Minkowski, Phys. Lett. B67, 421 (1977). 

41. J. Schechter and J. Valle, Phys. Rev. D25, 2951 (1982). 

42. The physics of Majorana neutrinos and 0vf3f3 are discussed in S. Bilenky and S. 
Petcov, Rev. Mod. Phys. 59, 671 (1987) [Erratum-ibid. 61, 169 (1987)]; 

B. Kayser, F. Gibrat-Debu, and F. Perrier, The Physics of Massive Neutrinos (World 
Scientific, Singapore, 1989); 

B. Kayser, Physica Scripta T121, 156 (2005). 

43. S. Pascoli and S.T. Petcov, Phys. Lett. B580, 280 (2003). 

44. Analyses of the possible values of | < turq > | have been given by H. Murayama and 

C. Pena-Garay, Phys. Rev. D69, 031301 (2004); 

S. Pascoli and S. Petcov, Phys. Lett. B544, 239 (2002); 

S. Bilenky, S. Pascoli, and S. Petcov, Phys. Rev. D64, 053010 (2001), and Phys. 
Rev. D64, 113003 (2001); 

H. Klapdor-Kleingrothaus, H. Pas, and A. Smirnov, Phys. Rev. D63, 073005 (2001); 

July 14, 2008 16:26 



26 1 . Neutrino mixing 

S. Bilenky et al, Phys. Lett. B465, 193 (1999); 
References in these papers. 

45. See also S. Elliott and P Vogel, Ann. Rev. Nucl. Part. Sci. 52, 115 (2002), and 
references therein. 

46. C. Albright, Phys. Lett. B599, 285 (2004). 

47. U. Seljak, A. Slosar, and P. McDonald, JCAP 0610, 014 (2006); 
J. Lesgourgues and S. Pastor, Phys. Rept. 429, 307 (2006). 

48. K. Mahn and M. Shaevitz, Int. J. Mod. Phys. A21, 3825 (2006). 

49. For an alternative approach to determining #atm 5 see M. Gonzalez- Garcia, M. Maltoni 
and A. Smirnov, Phys. Rev. D 70, 093005 (2004), and references therein. 

50. M. Fukugita and T. Yanagida, Phys. Lett. B174, 45 (1986). 

51. G. 't Hooft, Phys. Rev. Lett. 37, 8 (1976); 

V. Kuzmin, V. Rubakov, and M. Shaposhnikov, Phys. Lett. 155B, 36 (1985). 

52. S. Pascoli, S. Petcov, and W. Rodejohann, Phys. Rev. D68, 093007 (2003); 
S. Davidson, S. Pascoli, and S. Petcov, private communications. 

53. A. Cervera et al, Nucl. Phys. B579, 17 (2000); 
M. Freund, Phys. Rev. D64, 053003 (2001). 

54. If the magnetic moments are large, they might possibly tell us whether neutrinos 
are Majorana or Dirac particles. See S. Davidson, M. Gorbahn, and A. Santamaria, 
Phys. Lett. B626, 151, (2005); 

N. Bell et a/., Phys. Rev. Lett. 95, 151802 (2005); 

N. Bell et a/., Phys. Lett. B642, 377 (2006). 



July 14, 2008 16:26 



